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ABSTRACT 

We use near-infrared (NIR) spectroscopy obtained with GNIRS on Gemini, NIRSPEC on KECK, 
and ISAAC on the VLT to study the rest-frame optical continua of three 'Distant Red Galaxies' 
(having J s — K s > 2.3) at z > 2. All three galaxy spectra show the Balmer/4000 A break in the 
rest-frame optical. The spectra allow us to determine spectroscopic redshifts from the continuum with 
an estimated accuracy Az/(1 + z) ~ 0.001 — 0.04. These redshifts agree well with the emission line 
redshifts for the 2 galaxies with Ha emission. This technique is particularly important for galaxies 
that are faint in the rest-frame UV, as they are underrepresented in high redshift samples selected in 
optical surveys and are too faint for optical spectroscopy. Furthermore, we use the break, continuum 
shape, and equivalent width of Ha together with evolutionary synthesis models to constrain the age, 
star formation timescale, dust content, stellar mass and star formation rate of the galaxies. Inclusion 
of the NIR spectra in the stellar population fits greatly reduces the range of possible solutions for 
stellar population properties. We find that the stellar populations differ greatly among the three 
galaxies, ranging from a young dusty starburst with a small break and strong emission lines to an 
evolved galaxy with a strong break and no detected line emission. The dusty starburst galaxy has 
an age of 0.3 Gyr and a stellar mass of 1 xlO 11 M©. The spectra of the two most evolved galaxies 
imply ages of 1.3-1.4 Gyr and stellar masses of 4x 10 11 M Q . The large breaks for the two most evolved 
galaxies indicate that active galactic nuclei (AGN) do not dominate the rest-frame optical continuum 
emission of these galaxies, while for the younger starburst a significant contribution from an AGN 
cannot be ruled out. The large range of properties seen in these galaxies strengthens our previous 
much more uncertain results from broadband photometry. Larger samples are required to determine 
the relative frequency of dusty starbursts and (nearly) passively evolving galaxies at z ~ 2.5. 
Subject headings: galaxies: high redshift — galaxies: evolution — infrared: galaxies 



1. INTRODUCTION 

Recent studies have demonstrated that galaxies at z > 
2 show a large range in their rest-frame optical colors. 
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Deep near-infrared (NIR) imaging has allowed the iden- 
tification of a new class of z > 2 galaxies, complemen- 
tary to the Lyman break galaxies ( LBGs. iSteidel et"al] 
Il996albf) found in optical surveys. iFranx et all <|2003f) 
introduced the simple color criterion J s — K s > 2.3, 
to efficiently and successfully se lect red z > 2 galaxies 
(Distant Red Galaxies, or DRG s: lForster Schreiber et alJ 
2004; van Dokkum et afll2004D . The red colors can be 
caused by evolved stellar populations, dust, or a com- 
bination of both. It is difficult to assess the origin of 
the red colors from optical-to- NIR photometry alone (see 
iForster Schreiber et al.ll2004j) . The extension of photo- 
metric studies to the rest-frame NIR wavelength re gime 
helps: using IRAC on Spitzer, lLabbe et all p005|) dis- 
tinguish the old from the dusty system much more reli- 
ably than was previously possible, and find that a sig- 
nificant part (30%) of the DRG sample is indeed best 
described by old and passively evolving stellar popula- 
tion models. This result i s supported by the Spitzer 24 
(Ltm imaging presented bv I Webb et alJ J2006), who find 
that 65% of the DRGs host dusty star-forming stellar 
populations. Similar fi ndin gs have also been re ported by 
iPanovich et~afl l|2006|) and lReddv et all l)2005l) based on 
optical to mid-infrared photometry of the GOODS-South 
and North fields 

Although our insight in the nature of z > 2 galaxies has 
significantly improved by the recent access to the rest- 
frame NIR and IR wavelength regime, broadband photo- 
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metric studies have their limitations. First, as spectral 
details get lost, the origin of the observed emission is un- 
certain. Emission lines may affect the broadband fluxes, 
and a possible contribution by an active galactic nucleus 
(AGN) is difficult to quantify. Second, the modeling re- 
sults suffer from degeneracies between age, dust and the 
star formation history, especially when no spectroscopic 
redshift is available. 

The large samples of confirmed high-redshift galaxies 
availa ble today (e.g. St eidel et alJ l2003t iVanzella et alJ 
|2005|) give the impression that obtaining spectroscopic 
redshifts is relatively easy. The most efficient and pop- 
ular technique in obtaining redshifts for large samples 
of high redshift galaxies is multi-object spectroscopy 
at optical wavelengths. However, as about 75% of 
2 < z < 3 galaxies with M > 10 11 M Q have R > 25 
(van Dokkum et"aTll2006j) . most are well beyond the lim- 
its of optical spectroscopy. Thus NIR spectroscopy is 
needed to obtain redshifts for these UV-faint galaxies. 

For UV-faint passively evolving galaxies, which lack 
emission lines in their rest-frame optical, we are limited 
to the rest-frame optical stellar continuum to measure a 
redshift. The most easily detectable feature in the spec- 
trum of evolved stars, and thus the best option to confirm 
redshifts of z > 2 evolved galaxies, is the Balmer/4000 
A break. The break is also a powerful diagnostic in stel- 
lar population studies (e.g.lHamiltonlll985t iBalogh et alJ 
1999; iKauffmann et alJl2003a|) 

For galaxies beyond z ~ 1.5 the break shifts into the 
NIR, where the combination of high sky background and 
strong atmospheric absorption bands complicates con- 
tinuum studies of galaxies. Recently new NIR spec- 
trographs on 8-10 meter class telescopes have improved 
our access to the NIR regime significantly. This is 
demonstrated by the fi rst possible detections of rest- 
frame optical breaks by Ivan Dokku m et alJ <|2004f) and 
ISinmson et all pfiol with KECK/NIRSPEC and SUB- 
ARU/CISCO respectively. Furthermore, the large in- 
stantaneous wavelengt h coverage offered by GN IRS on 
Gemini-South (see e.g. Ivan Dokkum et aLll2005l) allows, 
for the first time, systematic continuum studies of galax- 
ies at z > 2. 

In this paper we present and study the rest-frame opti- 
cal continua of three DRGs, of which only one had a spec- 
troscopic redshift prior to the observations. The data 
are presented in §2. In §3 we fit the spectra with evolu- 
tionary synthesis models and compare the spectra with 
the broadband SEDs. In §4 we use spectral diagnostics 
to obtain stellar population properties. We present a di- 
rect comparison with other high-redshift and low-rcdshift 
galaxies in §5. We end with a summary and conclusions 
in §6. Throughout the paper we assume a ACDM cos- 
mology with Q m = 0.3, Q\ — 0.7, and H = 70 km 
s _1 Mpc -1 . All broadband magnitudes are given in the 
Vega-based photometric system. 

2. DATA 

2.1. Target selection and photometry 

The three studied galaxies lie in the field of MS1054- 
03, the CDF-South and the extended HDF-South re- 
spectively All fields have deep NIR imaging and 
accurate photometry in several optical-to-NIR bands. 
The MS1054-03 field is observed as part of the "ultra- 



deep" Faint InfraRed Extragalac tic Survey (FIRES) and 
the ph otometry is described bv iForster Schreiber et alJ 
(2006). T he Great Observatorie s Origins Deep Survey 
(GOODS; IGiavalisco et all l2004j) provides deep data of 
the CDF-South. The optical-to-NIR photometry that we 
used as part of this work will be described by Wuyts et 
al. (2006, in prep.). The photometry of the extended 
HDF-South is part of th e new MUlti-wavelen gth Survey 
by Yale-Chile (MUSYC; lOawiser et aT1l2005l Quadri et 
al. in prep). 

The three galaxies are chosen on the basi s of their K 
magn itude, red J s —K s color (J s — K s > 2.3: lFranx et alJ 
2003), and redshift between 2.1 and 2.7. To avoid a bias 
towards galaxies with strong UV emission, we did not re- 
quire that a previously measured spectroscopic redshift 
was available, and two out of three galaxies presented 
in this paper were selected on basis of their photomet- 
ric redshift. MS1054-1319 was the only galaxy which 
had a spectroscopic redshift (z — 2.423) prior to these 
NIR observati ons, derived from b oth t he Ly q and Ha 
emission lines l|van Dokkum et al.ll2003l 12004) . The NIR 
spectrum of CDFS -695 studied in this paper has already 
been presented by Ivan Dokkum et all ((21)05) , who give 
an emission line redshift of 2.225. 

Where possible the broadband fluxes were corrected for 
emission line fluxes. The observed infrared photometry 
for CDFS-69 5 is adjusted using t he lin e measurements 
presented by Ivan Dok kum et all p005|) . The emission 
line corrections are crucial to interpret the broadband 
photometry for this galaxy, as they account for 0.05, 0.17 
and 0.25 mag of the total J, H and K broadband mag- 
nitudes respectively. Hence, after corrections CDFS-695 
does not satisfy the J s — K s > 2.3 criterion. For this 
galaxy we have no observed optical spectrum. Fortu- 
nately, CDFS-695 is also part of the K20 survey, and 
from the spectrum presented bv lDaddi et al.l l|2004aD we 
conclude that the emission lines barely contribute to the 
observed optical broadband fluxes. Thus, no corrections 
are made to the broadband observed optical magnitudes 
of CDFS-695, but the photometric errors are increased by 
0.05 mag. For MS1054-1319 the emission line correction 
were derived bv lvan Dokkum et alJ (|2004) , who find con- 
tributions of 0.1, 0.02, 0.03, 0.04 mag to the broadband 
fluxes of B, J, H and K respectively, and an increase 
of the photometric error of V and I by 0.05 mag. For 
HDFS-5710, no emission lines are detected in the NIR 
spectrum and therefore no correction is applied. 

2.2. Observations 

Table 1 summarizes the observations. Two out of 
three galaxies presented in this paper were observed with 
GNIRS on Gemini-South in September 2004: CDFS-695 
and HDFS-5710 (program GS-2004B-Q-38). We used 
GNIRS in cross-dispersed mode, in combination with 
the short wavelength camera with the 32 1/mm grating 
(R=1000) and the 0'.'675 by 6'.'2 slit. In this configura- 
tion we obtained a wavelength coverage of 1.0 - 2.4 /xm, 
divided over 6 orders. Conditions on September 2 and 
3, during which we observed CDFS-695, were relatively 
good with a seeing of 077 in K and mostly clear sky. 
HDFS-5710 was observed on September 6 through inter- 
mittent clouds, and the seeing varied from 0'.'7 - l'/3. 

We observed the galaxies following an ABA'B' on- 
source dither pattern. In this way we can use the av- 
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TABLE 1 

Observations 



id 


K a 


V J s 


-K a 


Instrument 


date 


Grange 


Exp(min) 


Rspcc 


Seeing(") 


MS1054-1319 


19.01 


24.21 


2.58 


KECK/NIRSPEC 


2003/01/22 


2.05-2.47 


90 


1600 


0'.'9 














1.28-1.56 


60 


1600 


0'.'9 












2004/02/11 


1.23-1.55 


105 


1600 


O'.'S 














1.50-1.74 


60 


1600 


O'.'S 










VLT/ISAAC 


2004/04/25 


1.82-2.50 


81 


600 


Of.' 7 














2.14-2.26 


165 


3900 


Of.' 7 


CDFS-695 


19.12 


23.71 


2.33 


GEMINI-S/GNIRS 


2004/09/02-03 


1.0-2.4 


92 


1000 


0'.'7 


HDFS-5710 


19.31 


25.53 


2.47 


GEMINI-S/GNIRS 


2004/09/06 


1.0-2.4 


210 


1000 


l'.'O 



erage of the previous and successive exposures as the sky 
frame. The offset between A and A' ensures that differ- 
ent parts of the detector are used for sky subtraction. To 
optimize the signal-to-noise ratio (S/N) of the extracted 
one-dimensional (ID) spectrum, the shift between A and 
A' should be larger than the area over which the spec- 
trum is extracted. Compared to the standard ABBA 
pattern, this method decreases the noise by about 13 % 
in the sky-subtracted frames. 

We complemented the GNIRS spectra with the spec- 
trum of MS1054-131 9 obtained with N1RSPEC on Keck 
l|McLean et al] 11951 . We observed MS1054-1319 dur- 
ing two runs using the N4 filter and medium dispersion 
mode (1.23-1.55 /Ltm, R=1600), which is especially useful 
to target the region around the optical break for z ~ 2.5 
galaxies. The first run in January 2003 was character- 
ize d by cloudy weather and a typical seeing of 0'.'9 in 
K ifvan Dokku m et alJ l2004). The conditions during the 
second run in February 2004 were somewhat better with 
a seeing of 0'.'8. To extend the wavelength coverage for 
this galaxy, we complement the N4 coverage by H and 
K band spectra obtained with NIRSPEC on KECK and 
ISAAC on the VLT. For details on these observations 
see Table □ The NIRSPEC and ISAAC spectra were 
taken at three dither positions (ABC), as the slit length 
of both instruments is longer than for GNIRS. Similar 
to the ABA'B' pattern we can use the average of the 
previous and following exposures as sky frame. 

All targets were acquired by blind offsets from nearby 
stars. With NIRSPEC the alignment was checked and 
corrected if needed before every individual 900 s expo- 
sure. Furthermore, at each dither position a separate 
spectrum was taken of the offset star to determine the 
expected position of the object spectrum on the detec- 
tor. As GNIRS and ISAAC have better pointing stability 
there was no need for re-acquisition after each individ- 
ual exposure, and acquisition checks were performed ap- 
proximately every hour. The total exposure times of all 
objects and modes are listed in Tabled 

After each observing sequence, we observed an AV0 
star near the target. These spectra are used in the re- 
duction to correct for detector response and atmospheric 
absorption. 

We note that other targets were also observed during 
these runs. However, the galaxies discussed here are the 
only ones for which we specifically attempted to measure 
the continuum rather than just emission lines. 

2.3. Reduction of GNIRS spectra 



The reduction of t he NIRSPEC spectr a of M S1054- 
1319 is described bv Ivan D okkum et alJ l|2004[) . The 
ISAAC spectra of MS1054-1319 were reduced in a simi- 
lar way as the NIRSPEC spectra. For the GNIRS data 
reduction of CDFS-695 and HDFS-5710 we developed a 
suite of custom scripts to perform the following steps. 
We started the reduction by removing a bias pattern. 
The pattern repeats itself every 8 columns, differently 
in each quadrant. It was modeled by fitting a constant 
along the vertical direction in regions where no spectrum 
is present, and removed over the whole detector. We also 
corrected for persistence from the acquisition image on 
the detector. 

Cosmic rays were removed as follows. First, we re- 
moved the sky from the science frames by subtracting 
the average of the previous and next exposure. For addi- 
tional sky-subtraction we have to straighten the orders, 
which follow curved paths on the detector. We deter- 
mined the position of the object spectrum in each row 
of each order by tracing the spectrum of a bright star. 
Next, we straightened the orders by integer pixel shifts 
(to retain the cosmic ray shapes), using the expected ob- 
ject positions. Because of the small slit length, there 
is not enough 'empty' sky left to remove the remaining 
sky by fitting lines along the spatial direction. However, 
the negative object spectrum should cancel its positive 
counterpart, and any residual sky was removed by requir- 
ing t hat the biweight mean ijBeers. Flvnn. fc Gebhardtl 
11990ft along the spatial direction is zero at every wave- 
length. Next, we run L.A.Cosmic l|van Dokkuml 20011 
on the sky-free frames to identify the cosmic rays. The 
resulting cosmic ray masks were convolved by a boxcar 
to exclude neighboring pixels. The masks of the differ- 
ent orders were transformed back to the original shape of 
the spectrum without interpolation and combined to one 
cosmic ray mask for each individual exposure. Through- 
out the following description, the cosmic ray mask will 
be transformed in the same way as the science frames. 

We continued with the original raw frames and per- 
formed the same steps as described above to remove 
the sky, only this time we allowed interpolation when 
straightening the orders. For each exposure and each 
order, a ID noise frame was made from the subtracted 
sky. This noise frame was used to identify dead and hot 
pixels, which were added to the mask. 

Next, the orders of each frame were wavelength cal- 
ibrated using arc lamp frames. The calibrations were 
checked and if necessary corrected using skylines. Cor- 
rections for the instrumental response function and at- 
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Fig. 1. — Extracted one-dimensional spectra (upper panels) and their binned version (lower panels) for CDFS-695 (left), MS1054-1319 
(middle) and HDFS-5710 (right). The presented original ID spectra (upper panels) are smoothed by a boxcar of 25 A and 13 A (observed 
frame) for GNIRS and NIRSPEC spectra respectively. Bad atmospheric regions or wavelengths that are heavily contaminated by skylines 
are indicated in light grey. Bins that include emission line fluxes are plotted in grey as well. The best fit template spectra (dark grey) are 
overplotted (see Table |5J. For each galaxy a break is detected between 3650 and 4000 A. 



mospheric absorption were applied by dividing by the 
response spectrum. This spectrum was created from the 
observed spectrum of an AV0 star, which was divided by 
the spectrum of Vega. Residuals from Balmer absorption 
features in the spectrum of the AVO star were removed 
by interpolation. The spectrum of the AVO star was re- 
duced in the same way as the science frames. 

The different exposures were combined for each or- 
der, excluding cosmic rays and outliers as identified in 
the masks. Finally the orders were combined, properly 
weighting overlapping regions using the response func- 
tion. In the entire procedure pixels were interpolated 
only once, to minimize smoothing and noise correlations 
in the final frames. 

2.4. Extraction of 1 dimensional spectra 

For all three galaxies the ID spectrum was extracted 
by summing all lines (along the wavelength direction) 
with a mean flux > 0.25 x the flux in the central row, 
using optimal weighting. In addition, a 'low-resolution' 
(binned) spectrum was extracted as follows. The pixels 
were sampled along the wavelength direction in each line 



of the two-dimensional (2D) spectrum using the biweight 
estimator. Thus we created a binned spectrum for each 
line of the 2D spectrum. Hereafter the lines with a mean 
flux > 0.25 x the flux in the central line were added 
corresponding to their weighting factors. This procedure 
gives a higher S/N in the final low resolution spectrum 
than binning the original extracted ID spectrum. 

Regions with low or variable atmospheric transmis- 
sion or with strong sky line emission were excluded. For 
CDFS-695 and HDFS-5710 observed with GNIRS, we de- 
fined these regions as those with less than 5% of the max- 
imum transmission and with sky line intensities of 30% 
or more of the strongest line. For MS1054-1319 observed 
with NIRSPEC and ISAAC, these criteria were <30% 
and >25%, respectively. For MS1054-1319 the spectra 
were sampled in 30 pixel bins (~80 A), and the other 
two galaxies were sampled in bins of 50 pixels (^250 A). 
Rather than to vary the number of pixels that contribute 
to each wavelength bin (i.e. excluding bad wavelength 
regions), the width was adjusted such that each bin con- 
tains 30 or 50 'good' pixels. 
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For MS1054-1319 we combined the NIRSPEC and 
ISAAC spectra of different but overlapping wavelength 
coverages. Adding 2D spectra would have decreased the 
S/N as there is a difference in seeing and sampling of 
instruments. Instead we combined the different spectra 
in the following way. First, we extracted a low resolu- 
tion spectrum for each individual 2D spectrum, by us- 
ing the same bins and a common mask for excluding 
regions of poor atmospheric transmission or with strong 
sky lines. These spectra were scaled by minimizing \ 2 m 
the overlapping regions and thereafter combined. These 
scaling factors were used to combine the 'high resolu- 
tion' ID spectra as well. Spectra of different wavelength 
coverages without overlapping regions were scaled using 
the total broadband H and K magnitudes. To scale the 
spectra, we modeled the continuum by fitting a straight 
line to the low resolution spectrum, excluding bins that 
may contain emission lines. The modeled continuum was 
converted to F v and thereafter convolved with the appro- 
priate filter-curve and normalized using the emission line 
corrected total fluxes. The GNIRS spectra of HDFS-5710 
and CDFS-695 were scaled in the same way, but as the 
wavelength coverage is not divided over different bands, 
we used an error-weighted average scale factor obtained 
from the broadband H and K magnitudes to tie the spec- 
trum to the broad band photometry. In this procedure, 
we did not use the broadband J magnitude because the 
Balmer/4000 A break complicates the modeling of the 
continuum in this band. 

The ID original and low resolution spectra are pre- 
sented in Fig.Q] For MS1054-1319 we only show the area 
around the break observed with NIRSPEC, which has 
continuous coverage. For each galaxy an optical break is 
detected. This is no coincidence, as the selection crite- 
rion of J s — K s > 2.3 was introduced to identify galaxies 
with prominent 4000 A and/or Balmer breaks. There 
are even indications that some individual absorption fea- 
tures are detected for MS1054-1319 (H( and Ca(H)) and 
CDFS-695 (Ca and H<5). The sensitivity of the spectra 
of CDFS-695 and HDFS-5710, observed with GNIRS, de- 
creases in the blue, which is reflected in the error bars 
in Fig. The outliers bluewards of the break in CDFS- 
695 are caused by residual bias patterns which cannot be 
properly removed. 

In the spectrum of CDFS-695, we clearly detect several 
emission lines that are labeled in Fig. [T] The analysis of 
these emission lines is presented by Ivan Dokkum et aTl 
( 2005). For MS1054-1319 we observe two possible emis- 
sion lines ([Oil] at 3727 A and [0 III] at 5007 A) in the 
combined NA and H band spectra. This galaxy also has 
detected Ha and [Nil] em ission in its K-band spectrum 
Ijvan Dokkum et al.ll2004[) . Strikingly, for HDFS-5710 no 
emission lines were detected at all. Note that the galaxy 
without emission line detections appears to have the 
strongest break (HDFS-5710), and that the galaxy with 
the strongest lines (CDFS-695) has the weakest break. In 
§ 14.21 we will investigate this relation between the break 
strength and emission lines in more detail. 

3. SPECTRAL MODELING 

Modeling the broadband photometry is a popu- 
lar method to study the properties of high-redshift 
galaxies (e.g. Sawicki k, Yedll998t iPapovich et aLll2001i 
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Fig. 2. — Broadband photometry together with the observed 
and modeled spectra for CDFS-695 (upper panel; B, V, I, z, J, H, 
& K), MS1054-1319 (middle panel; U, B, V, R, I, J, H, & K), 
and HDFS-5710 (lower panel; U, B, V, R, I, z, J, H, & K). The 
error bars represent the lcr uncertainties of the flux measurements. 
The bins that are contaminated by emission line fluxes, and thus 
not included during fitting are plotted in grey. The best fits to 
just the photometry (light grey) and to the spectra in combination 
with the rest-frame UV bands (dark grey) are overplotted. For 
MS1054-1319 and CDFS-695 the redshift was fixed to the emission 
line redshift during the fitting procedure, while for HDFS-5710 z 
was a free parameter. The r, A v , and age (see §3.1) of the best fit 
are listed in the panels. 



fortunately, the modeling results suffer from degeneracies 
between age, dust and the star formation timescale. In 
this section we investigate the additional constraints pro- 
vided by the rest-frame optical spectra. 



iShaplev et al.ll2001t iForster Schreiber et al.ll20oH) . Un- 



3.1. Fitting procedure and results 

used synthetic spectra by iBruzual fc CharloU 
( 2003) to m odel our data. We selected a stellar popu- 
lation with a lSakoeterl l)1955|) initial mass function (IMF) 
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TABLE 2 
Modeling Results 



id 


T 

Gyr 


age 
Gyr 


A v 
mag 




2 b 


xL 


xL N c 


10 11 M 


SFR 

Moyr" 1 


SFR/Af* 
10~ n yr" 1 


CDFS-695 
phot 


03+ 9 ' 97 


o 05+ 015 
u - uo -o.oo 


1 q+0.2 




0.77 


1.46 


3.46 3 


0.8l° j 


799 +874 
'^-662 


Q o 9 +1276 
ao -854 


spec 


o 25+ 9 - 75 

U ' ZO -0.22 


oq+0- 28 


1 4+ 01 




2.78 


2.85 


2.95 39 


1 n + 0.3 

L -°-0.3 


9S4+ 100 
2»4_ 206 


226:% 



MS1054-1319 

phot 0.50^'jg l-6lto.I? 0.7tg'.? - 1-28 2.10 3.20 4 4.2±]J 47±J 1 111? 

spec 0.40l°J° l-28l°;^ 0.7±|^ - 1.59 1.67 1.73 61 3.5±g-_| 48±^ U±% 

HDFS-5710 

phot O.OSt^s 0.20^;?° 1.9^ 2-28±°:|g 0.52 1.90 2.83 4 3.2±»;J 4±* 417 

spec 0.32±g-^| 1.43±J-|| 0.9±J'g 2.283±gg|3 2.21 2.37 2.57 38 4.0^1 ^t\f 4±| 3 

a The redshifts of CDFS-695 and MS1054-1319 arc fixed to z Ha =2.225 and z H c=2.423 respectively. 

^X 2 i s given per degree of freedom. The la and 2a values correspond to the 68% and 95% confidence levels, derived from the Monte 
Carlo simulations described in §3.1 
c Numbcr of degrees of freedom 



between 0.1-100 M Q , a solar metallicity, and based on 
the Padova 1994 evolutionary t racks . We adopted the 
reddening law of iCalzetti et alJ l)2000f) and corrected for 
intergalactic attenuation by the Lyq forest using the pre- 
scriptions of lMadau et alJ l)1996f ). For a discussion about 
the choice of t he model parameters and the eff ects of their 
variations, see lForster Schreiber et all l)2004[) . 

The star formation history is parametrized by an ex- 
ponentially declining function with characteristic time 
scale r. Synthetic spectra were constructed for a grid of 
31 different r's between 10 Myr and 10 Gyr, 24 different 
ages between 10 Myr and 3 Gyr (and always restricted to 
be less than the age of the universe at a given redshift), 
redshifts in steps of Az = 0.001, and 31 extinction val- 
ues (A v ) between and 3 mag. For MS1054-1319 and 
CDFS-695, the redshift was fixed to the emission line 
redshift. 

The templates were sampled in the same bins as the 
observed spectra, thus wavelength regions with bad at- 
mospheric transmission and strong skylines in the ob- 
served spectra were excluded in the synthetic spectra as 
well. Additionally, bins that were contaminated by emis- 
sion line fluxes ([OIII], Ha, [Nil] and [SII] ) were excluded 
during fitting, as the models do not include emission by 
the ionized interstellar medium. To extend our wave- 
length coverage to the rest-frame UV, we complemented 
the observed spectra with the optical photometry. The 
broadband fluxes of the synthetic spectra were derived by 
integrating the flux accounting for the filter curves. The 
minimum error in the broadband fluxes was set to 0.05 
mag to account for absolute calibration uncertainties. 

The spectra in combination with the rest-frame UV 
photometry were fitted by using least-square minimiza- 
tion (x 2 )- We performed 200 Monte Carlo simulations to 
calibrate the confidence levels. For these simulations we 
varied the fluxes randomly within the errors, assuming a 
Gaussian distribution, and followed the same procedure 
as applied to the original data. The la and 2cr confidence 
levels are derived from the x 2 values of the original fit 



whic h enclose the best 68% and 95% of the simulations 
(see iPanovich et alJl2"00lD . These x 2 values per degree 
of freedom are listed in Table 2. 

The best spectral fits are shown in Fig. [21 The r, age, 
Ay, redshift, stellar mass, star formation rate (SFR) and 
X 2 corresponding to these fits are listed in Table El All 
allowed combinations of z, t, Ay and age are presented 
as la and 2cr confidence levels (color contours) in Fig.|3 

3.2. Continuum redshifts 

To test the accuracy of the redshift obtained by spec- 
tral fitting, both CDFS-695 and MS1054-1319 were re- 
fit with z as a free parameter. The spectra of these 
galaxies give remarkably good solutions for z. This 
is illustrated in Fig. H For MS1054-1319, which has 
an Ha and Lyq redshift of 2 . 423 a nd 2.424 respec- 
tively l|van Dokkum et all 120031 12004th we find z fit = 
2.423±£;E$ 2 . For CDFS-695 w hich has an Ha and Lyq 
redshift of 2.225 and 2.223 ijvan Dokkum et alJ 120051 
iDaddi et aH2004a|) . we find z& t = 2.21±g-J|. This galaxy 
has a secondary solution for redshift (z ~ 2.32) that is 
probably caused by mis-identification of the absorption 
features in this spectrum. 

The uncertainty on the derived value for the redshift is 
decreased significantly compared to only using the broad- 
band fluxes. The green lines and shaded areas in Fig. 0] 
indicate t he la redshift range allow ed by the photome- 
try, using iBruzual fc Charlotl (|2003) models (see confi- 
dence levels in Fig. 0). The purple lines and surround- 
ing shaded areas are obtained by fitting linear combi- 
nations of theoretical and empirical templates to the 
br oadband photometry, follow ing the method described 
hv lRudnick et alJ (j2001i 12003ft . The uncertainty on the 
redshifts are decreased by a factor of 6-50 when including 
the spectrum. 

For HDFS-5710 which is the only galaxy without an 
emission line redshift, spectral modeling yields z&t = 
2.283±g;g|g (see Fig. gj). This result illustrates the po- 
tential of NIR spectroscopy to measure redshifts of z > 2 
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Fig. 3. — Comparison of modeling results between the spectrum plus rest-frame UV broadband fluxes and the broadband SED for CDFS- 
695 (upper panels), MS1054-1319 (middle panels) and HDFS-5710 (lower panels). The best-fit solutions and their lc and 2cr confidence 
levels are presented by the crosses and color contours for the spectrum together with the optical photometry, and pluses and filled greyscale 
contours for solely the optical-to-NIR photometry. For MS1054-1319 and CDFS-695, the redshift was fixed to zjja to derive r, age, Ay, 
SFR and stellar mass. The dotted lines in the age versus z plots indicate the emission line redshifts. For HDFS-5710 the redshift was a 
free parameter when deriving all properties. The \ 2 value of the best fit solution is given in the left panels, in red for the spectral fit and 
in black for the photometry fit. This figure shows that the spectrum improves the constraints on all properties substantially. In addition, 
the modeled redshifts for CDFS-695 and MS1054-1319 are remarkably accurate. 



galaxies that are faint in the rest-frame UV and have no 
detectable emission lines. 

3.3. Comparison to broadband SEDs 

To investigate if the spectra confirm and improve the 
constraints on the stellar populations properties, we re- 
fitted the photometry of the three galaxies, ignoring the 
spectral information. The comparison between the al- 
lowed solutions for the input parameters z,T,age and 
Ay , and the output parameters stellar mass and SFR for 
the spectra and broadband SED is presented in Fig. [3J 
The corresponding values are listed in Tabled The best- 
fit templates to the spectra and photometry, and the best 
fits to the photometry alone are shown in Fig. [5] 

Fig-EDshows that the spectra greatly improve the con- 
straints on the modeled properties, for all three galaxies. 
Most of the la solutions to the spectra fall within the 
2a confidence levels of the photometry. However, most 
of the la confidence regions of the photometric fits are, 
within 2a, not allowed by the spectral fits. For example, 
while the photometry of CDFS-695 still allows a wide 
range in dust content and age, the solutions dustier than 
1.6 mag and younger than 0.1 Gyr are ruled out by the 
spectrum. Even MS1054-1319, for which the properties 
are already well-constrained by the photometry, the spec- 
trum decreases the la confidence interval of age and r 
by a factor 7-8. And finally, for HDFS-5710 the largest 
gain is achieved by the better constrained redshift. Also, 



solutions with r longer than 700 Myr are ruled out by 
the spectrum. 

The normalized x 2 values per degree of freedom for the 
fits to the broadband photometry are lower than those 
for the fits including the NIR spectroscopy. However, 
the difference between the x 2 of the best fit and the la 
and 2a confidence intervals are smaller for fits includ- 
ing the NIR spectra (see Table |3J), which result in the 
tighter constraints. One possible cause for the difference 
in minimal x 2 of the best photometric and spectral fits 
is that x 2 statistics are very sensitive to the data uncer- 
tainties. The photometric errors may be overestimated 
or the spectral errors may be underestimated. Another 
cause could be template mismatch, which can become 
more problematic with more detailed spectroscopic in- 
formation. 

To test the robustness of our results we repeated the 
fitting procedure for MS1054-1319, artificially decreas- 
ing the errors on the photometry by a factor of 2 and 
increasing the errors on the spectra by a factor of 2. The 
allowed ranges for Ay, age and r are slightly increased 
and reduced for the spectra and photometry respectively, 
but even in this extreme case inclusion of the spectra still 
improves the constraints on the stellar populations. 

3.4. Discussion of modeling results 

Our fitting results show that including the NIR spectra 
when fitting the broadband fluxes sets tighter constraints 
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Fig. 4. — Comparison between continuum redshift, broadband 
photometric redshift and the emission line redshift. The best-fit 
redshift and its la confidence interval to the NIR continuum to- 
gether with the optical photometry is indicated by the red line and 
surrounding shaded area respectively. The green thick line and 
the surrounding shaded area present the best-fit redshift and its 
allowed la range when fitting the op tical-to-NIR broadband pho- 
tometry by Bruzual & Chariot (2003) models. The best-fit redshift 
and its la uncertainty obtained when fitting the photometry by 
linear combinations of emp irical and theoretical galaxy templates 
(Rudnick et al. 2001, 2003) is shown by the purple thick lines and 
shaded areas. The emission line redshift, if present, is indicated 
by the black dotted line. This diagram illustrates that the contin- 
uum shape is powerful in obtaining accurate and well-constrained 
redshifts. 



on stellar population properties and reduces degeneracies 
between age, extinction and the SFH. This is mainly due 
to the measurement of the shape and strength of the op- 
tical break. We conclude that the spectral shape in com- 
bination with the photometry is more powerful than the 
broadband photometry alone in tracing the properties of 
stellar populations in high redshift galaxies. 

The DRGs studied here span a wide range in their 
derived properties, consistent with results of larger sam- 
ples based on modeling of th eir broadband SEDs (e.g. 
IForster Scnreiher et al.ll2004t ILahhe et al .1120051). Their 
SFH, reddening and age are quite different, ranging from 
CDFS-695 which is about 0.2 Gyr old and is forming 
a few hundred solar masses a year, to HDFS-5710 and 
MS1054-1319 which arc both best fitted by a stellar pop- 
ulation of 1.4 Gyr old and are forming a few tens of solar 
masses a year. 

lLabbe et alJ i|2005D divided their DRG sample into 
dusty starburst and "red and dead" galaxies based on 
the observed NIR and IR colors (I — K and K — 4.5/i). 
T o compare the DRG s presented in this work with those 
of [Labbc eTal] l|200l . we measured the same colors by 
extrapolating the best spectral fit to the rest-frame NIR. 
All three galaxies fall in the dusty star-forming part of 
the plot, although their locations move closer towards 
the "red and dead" part, when going from CDFS-695 to 
MS1054-1319 to HDFS-5710. 

4. SPECTRAL DIAGNOSTICS 

In the previous section we studied the spectra by com- 
paring them to models. The break appeared to be 
the driving feature behind the improved modeling con- 



straints. In this section, we take a different approach 
by measuring the break directly, and use this feature in 
combination with Ha to constrain the star formation his- 
tories of the galaxies. 

4.1. Break indices 

The Balmer and 4000 A break are often treated as one 
feature, owing to their similar locations and the fact that 
they partially overlap. However, the breaks originate 
from different physical processes, and behave differently 
as populations age. Both breaks are due to absorption in 
the atmosphere of stars. The 4000 A break arises because 
of an accumulation of absorption lines of mainly ionized 
metals. As the opacity increases with decreasing stel- 
lar temperature, the 4000 A break gets larger with older 
ages, and is largest for old and metal-rich stellar popula- 
tions. The metallicity is of minor influ ence for ages less 
than 1 Gyr ijBruzual k, Char lotj 12003]) . There are sev- 
eral definitions introduced to q uantify the strength of the 
4000 A break. iBruzuall l|1983H proposed L>(4000), which 
measures the ratio of the average flux density F v in the 
bands 4050-4250 A and 3750-3950 A around the break. 
Because of the broad regions, this index is f airly sensitive 
to red dening by dust. To reduce this effect . IrTalogh et ail 
(1999) defined a new index D„(4000), based on smaller 
continuum regions (3850-3950 A, 4000-4100 A). 

The Balmer break at 3646 A marks the termination of 
the Hydrogen Balmer series, and is strongest in A-type 
stars. Therefore the break strength does not monotoni- 
cally increase with age but reaches a maximum in stel- 
lar populations of intermediate ages (0.3-1 Gyr). The 
strength of the Balmer sequence can be best measured 
from the individual Balmer lines, such as H5. However, 
as our spectra do not allow the measurement of this fea- 
ture, we use the strength of the Balmer break (Db), 
which we define as the ratio of the average flux density 
F v in the bands 3500-3650 A and 3800-3950 A around 
the break. The large regions are not optimal, but are a 
trade-off between dust-dependence and having sufficient 
S/N using the spectra of high redshift galaxies. This 
index is also partially influenced by the 4000 A break. 
Nevertheless, the age- dependence of Db is very similar 
to that of Hi5. 

We determined D n (4000) and D B for the three DRGs. 
The average flux in the regions around the breaks is 
measured in the same way as we extracted the binned 
spectra. For MS1054-1319 it was impossible to mea- 
sure the flux redwards of the 4000 A break, as the 
defined region is entirely covered by an atmospheric 
band. Therefore we me asured 13(4000) and used the 
iBruzual fc Chariot! l|2003fl models to derive £>„(4000) for 
this galaxy. When converting the indices we allowed an 
Ay between and 3 and all three star formation his- 
tories to determine the uncertainty. We note that this 
only slightly increases the uncertainty as for the mea- 
sured break strength of MS1054-1319 the correlation be- 
tween 13(4000) and Z3„(4000) is not very dependent on 
dust or the star formation history. 

The break measurements are listed in Table |21 and 
shown in Fig. The behavior of the Balmer and 4000 
A break with age for TioMyn 7"320Myr and TioGyr (which 
is comparable to constant star formation for galaxies at 
z > 2) models are overplotted for the first 3 Gyr. This 
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TABLE 3 

Spectral features 



id 


W HQ (A) a 


D„(4000) 


D B 


CDFS-695 


99±10 b 


i 17+0.14 


1 84+ ' 24 


MS1054-1319 


19±4 C 


i 9H +0.16 


1 69+ 024 
l.o»_ 2 Q 


HDFS-5710 


< 10 


i ,0+0.47 
L - oa -0.33 


2 40+ 67 
^■* u _0 .43 



a Uncorrected for Balmer absorption 
trn 



var^Dokkvrrri^^alJ 120^ 
wtn^DcTtkumet^uJ 12004 



figure is useful for discriminating between different star 
formation histories, as it provides a clear illustration that 
(near-)constant star formation models are not able to 
produce a large break within 3 Gyr. However, if a galaxy 
is extremely dusty, reddening can mimic a very large 
break. All three galaxies fall on the model curves. For 
MS1054-1319 the break implies an age of about 1 Gyr. 
The break of CDFS-695 reveals a younger stellar popu- 
lation. The large Balmer break for HDFS-5710 suggests 
that this galaxy is in a post-starburst phase. Unfortu- 
nately the errors are too high to draw firm conclusions 
from these diagnostics. 

4.2. Comparison to Ha equivalent width 

In addition to the break, the equivalent width of Ha 
(Who) can help discriminate between dusty starbursts 
and evolved stellar populations. As Wh q is sensitive to 
the ratio of the current and past star formation it is an 
independent measure of the star formation history and 
the age of a stellar population. 

The W h» of MS1054-1319 and CDFS- 695 were mea- 
sured bv Ivan Dokkum et all l)2004 l2005[h For HDFS- 
5710, no Ha emission line is detected, and thus we give 
an 3cr upper limit for Wu a using the redshift derived in 
Sect. El and assuming a FWHM of 500 km/s. To check 
for a possible Ha detection within the 2a allowed red- 
shift range, we measured the S/N of a possible line for 
redshifts between 2.2 and 2.5 (in steps of 0.001). For 
none of these redshifts do we find a detection exceeding 
the limits derived for the best fitting redshift. The values 
of Whq for all galaxies are listed in Table |2| 

In Fig.|HJi we plot Wh« as a function of £>„(4000) for 
the galaxies and different models (rsoMyr, 7"320Myr and 
TioGyr)- The plotted Whq are corrected for a Balmer 
absorption with an equivalent width of 4 ± lA, which 
is typical for stellar populatio ns with ages of 0.5-2.0 Gyr 
ijForster Schreibe r et al.l 2004D. T he model tracks arc de- 
rived from the iKerrnlcuttl l)1998j) law, which relates the 
luminosity of Ha t o the S FR, in combination with the 
iBruzual fc Chariot! l)2003j) models. These tracks illus- 
trate that neither a large break nor a low Wh q can be 
produced for Tiogyr or CSF models within 3 Gyr. 

Fig. illustrates that the three galaxies are at dif- 
ferent stages in their stellar evolution. The undetected 
Ha of HDFS-5710 is consistent with the large break, as 
both imply an evolved stellar population. Furthermore 
we note that the SFR from the continuum fits (see t |3.1fl 
of 19 Mo/yr is in agreement with the SFR upper limit 
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Fig. 5. — The Balmer break versus the 4000 A break for the 
three studied galaxies and three different star-formation models. 
The small black diamonds and the corresponding values indicate 
the age in Gyr. The effect of attenuation is indicated as a vector in 
the plot. Measurem ents for three type of AGNs are shown as well 
IFrancis et all 199 II) . This figure illustrates that a dust-free rioGyr 
(which is comparable to constant star formation for galaxies at 
z > 2) is unable to produce large breaks within 3 Gyr. For MS 1054- 
1319 the combined break implies a stellar population of about 1 
Gyr, while for CDFS-695 we find a younger stellar population. The 
large Balmer break of HDFS-5710 suggests that this galaxy is in 
a post-starburst phase, and that riogyr or CSF is highly unlikely. 
The large breaks for MS1054-1319 and in particular HDFS-5710 
imply that the optical spectrum is dominated by stellar light. 



of 25 MQ/yr, derived from the 3er upper limit on Ha and 
assuming the best-fit Ay of 0.9 mag. For MS1054-1319 
Who provides additional evidence that the stellar light 
is dominated by evolved stars. And, in agreement with 
all the modeling results, CDFS-695 is still forming stars 
at a high rate. 

We note that for CDFS-695 and MS1054-1319 there 
are indications that Ha is not entirely due to photo- 
ionization from star formation. The unusually high 
[Nil] /Ha ratios could suggest t hat another source of ion- 
izatio n is contributing to Wh l)van Dokkum et all200i 
2005) of these two galaxies. 

4.3. Active Galactic Nuclei? 

Optical spectroscopy, X-ray studies, and 24/im imag- 
ing infer that about 5% to 25% of the DRGs show signs of 
AGN activity (Wuyts et al. in prep . ; iRubin et alJl200l 
Reddv et all 120051 iPapovich et"all l2006t iWebb et alJ 
20061) . depending on if -band depth of the sample and 
the type of AGN. AGNs affect rest-frame optical spec- 
tra in two ways; they contaminate the stellar continuum 
emission and produce characteristic emission line shapes 
and ratios. Continuum contribution from an AGN weak- 
ens the stellar break strength, and the strong breaks for 
HDFS-5710 and MS1054-1319 (see Fig.0 imply that an 
AGN cannot be the dominant contributer to the rest- 
frame optical continua of these galaxies. For CDFS-695 
the observed optical break is relatively weak and a sig- 
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Fig. 6. — Comparison of spectral properties of the thre e studied galaxies with LBGs (green, Erb ct al. 2003), and SDSS galaxies (small 
grey dots, Tremonti ct al. 2004; Kauffrnann et al.|2003a|bl). Evolution tra cks for tioq vi . (so lid line), T32QMyr(dotted line) and t.-,om v , (dashed 
line) models, derived from the Bruzual & Chariot (2003) library and the Kcnnicutt (1998) law are overplotted without correction for dust. 
The small black diamonds and the corresponding values indicate the age in Gyr. The correction for the integrated attenuation of the whole 
galaxy and extra extinction towards HII regions on the models are indicated by vectors, (a) Wjja versus D n (4000). (b) Whq versus the 
stellar mass, (c) D n (4000) versus stellar mass. The masses of the galaxies are derived from the best-fit stellar model. 



nificant AGN contribution to the continuum cannot be 
ruled out. We note that a minor nuclear contribution, 
as for a Seyfert 2 (Fig. [SJ , would imply that the stellar 
breaks are even stronger, which would yield older ages 
for the galaxies. 

There is additional evidence that MS1054-1319 and 
CDFS-695 do not host an AGN which dominates the 
rest-frame optical continuum. NICMOS imaging of 
MS1054-1319 shows that the emission from this galaxy 
is extended (Toft et al. in prep). Furthermore, this 
galaxy is not detected in X-rays ( Ru bin et alJ I2004j) 
and there is no in dication for an AGN in the rest- 
frame UV-spectrum (jvan Dokkum et aTllkXM ID 1671). 
For CDF S-695 the res t -frame UV spectrum lacks AGN 
features (IDaddi et al.l I2004a|) . although the unusual 
high [Nil] /Ha ratios mentioned in the previous sec- 
tion could indicate an AG N contribution to this galaxy 
Ijvan Dokkum et alJl2005^ . 

5. COMPARISON TO OTHER GALAXIES 

In the previous sections we learned through modeling 
and measuring the spectra that the DRG sample contains 
galaxies with a wide variety of properties. In this section 
we compare the measured and modeled properties with 
those of other high-z and low-z galaxy samples. 

5.1. Comparison to other high-redshift galaxies 

Several studies find that at similar rest-frame V 
magnitude and redshifts, DR Gs are older, more mas- 
sive, and dustier than LBGs (van Dokkum et alll200l 
iForster Schreiber et all l200l ILabbe et al.M2005|) . We 
compare in this section the properties of the studied 
DRGs with those of LBGs at similar redshift, using the 
spectral indices and modeling results derived in previ- 
ous sections. Ideally the samples should be matched in 
K band magnitudes or stellar masses. As such matched 
samples are not available caution is required when inter- 
preting the differences. 

As there are no optical break measurements of LBGs, 
we have to rely on the 'predicted' break, based on the 



best fitted SEDs. We used t he BX and MD galaxies pre- 
sented by lErb et all (J2003) , which have a median red- 
shift of z ~ 2.3. The photometry and best fit models of 
the galaxies are kindly provi ded by D.K. Erb; 4 galax- 
ies of these are presented bv lShaplev et alJ ((2005) . The 
break has been derived for each of the 12 galaxies in the 
Q1623 and Q1700 fields individually and yield a median 
D n (4000) of 1.21+o;o6. The W Hq for th e BX/MD s ample 
is derived from the Lhq presented bv lErb et al.l (2003) 
and the best fit models. These values are quite uncertain 
as slit-losses can int roduce errors of up to about 50% in 
the absolute fluxes l)Pettini et al.ll200l|) . The predicted 
D n (4000) and W Ha of the BX/MD sample are plotted 
in Fig. Hju 

While our galaxies span different stages of stellar pop- 
ulation evolution, the LBGs are located on the young 
and star forming part of the plot (Fig. [Hla.) . Both DRGs 
and LBGs fall nicely on the predicted T320Myr track. 

The LBGs plotted here are all fainter than K = 20, and 
thus it could be that brighter LBGs ha ve properties more 
similar to those of the pl otted DRGs llAdelberger et ID 
l2005HShaplev et al.ll2?)f)l . Also the difference in stellar 
mass of the LBGs and DRGs could be the cause of the 
different properties found for both samples. To show how 
both properties relate to stellar mass, we plotted Wh q 
and -D„(4000) versus the stellar mass in Fig. [fj) and|SJ: 
respectively. The DRGs MS1054-1319 and HDFS-5710 
do not only have larger br eaks and smaller Wh« than 
the BX/MD galaxies from lErb et all (|2003j) . but have 
also much higher stellar masses. As they span different 
ranges in stellar mass, it is difficult to directly compare 
the properties of the published samples, or interpret it in 
the context of evolutionary scenarios linking LBGs and 
DRGs. 

5.2. Comparison with low-z galaxies 

To examine the implications for the evolutionary his- 
tories of nearby galaxies we now compare the breaks, 
Whu, and masses of z > 2 galaxy populations with 
those of SDSS galaxies. The L>„(4000) and W Hct (uncor- 
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rected for dust) of the S PSS galaxies are adopted from 
the catalogs pr e sented in iKauffmann et al.1 l)2003al) and 
iTremonti et alJ l)2004|) respecti vely. We used the media n 
dust-corrected stellar ma sses (jKauffmann et all l2003aj) . 
Kauffmann_el_alJ l2 QQ3aD used the same stellar library 
(Bruzual & Chariot! l2003[) to deriv e the stellar masses, 
but adopted the IMF by iKroupal ll2001fl At a given 
rest-frame V-band luminosity, the IKroupal 1)200 ID IMF 
yields masses which are a bout a factor of 2 lower than 
the masses o btained using a iSalpeterl lll955|) IMF between 
0.1-100 M Q ijBruzual fc Charlotll2003|) . and we applied 
this correction to the SDSS galaxies. 

In Fig. the SDSS galaxies show a tight relation be- 
tween D„(4000) and W Hq - Remarkably, DRGs and LBG 
fall on this same low-redshift relation. The similarity of 
the z ~ 2.3 galaxies and SDSS galaxies breaks down 
when we include the stellar mass in Figs. EJd and c. At a 
given mass the z ~ 2.3 galaxies have on average smaller 
breaks and higher Whc« than the low-redshift galaxies. 
This implies that they are younger and have higher spe- 
cific SFRs than their low-redshift analogs. Furthermore, 
galaxies with the properties observed for the z ~ 2.3 
DRGs are rare at low redshift. 

To show how simplified SFHs for HDFS-5710 and 
MS1054-1319 behave in this plot, we plotted the tracks 
for the T320Myr and TioGyr models, scaled such that both 
models form the stellar mass of MS1054-1319 and HDFS- 
5710 within 3 Gyrs. A straightforward explanation for 
the location of DRGs on Fig.EK~c is that they evolve into 
galaxies with properties comparable to those of the most 
massive and oldest galaxies in the SDSS. This may imply 
that the DRGs are the younger progenitors of galaxies of 
similar mass today. 

6. SUMMARY AND CONCLUSIONS 

We have presented rest-frame optical spectra of three 
DRGs at z ~ 2.3. To avoid a bias towards galaxies 
with strong UV emission, we selected objects which did 
not necessarily have a spectroscopic redshift from emis- 
sion lines prior to the observations. We were able to 
measure the continuum shape and directly detect the 
Balmer/4000 A break for all galaxies we observed for 
this purpose. This presents a significant advance be- 
yond previous studies whic h were based on broadband 
photometry (e.g. I Papovic h et alJ 12 001: Shaolcv e£jiL| 
[ 20011 iRudnick et alJl2003t iForster Schreiber et alJl2004i 
IPapovich et all 120061) or r e st-frame optical emission 
lines (e.g. lErb et all 12003 Ivan Dokkum et all 12004 
Swin bank et alJl2004|) . 

We have explored how the direct measurement of the 
optical continuum shapes and in particular the Balmcr 
and 4000 A break can contribute to our understanding 
of z > 2 galaxies. First, we have demonstrated that 
fairly accurate redshifts for red z ~ 2.3 galaxies can 
be obtained by modeling the rest-frame optical contin- 
uum around the Balmer/4000 A break. The uncertain- 
ties on the modeled redshifts range from Az/(1 + z) ~ 
0.001 — 0.04, depending on the S/N of the spectrum and 
the strength of the break. Our continuum redshifts agree 
very well with the emission line redshifts available for 2 
of the 3 galaxies (within 0.5%). This opens up the possi- 
bility of deriving redshifts for UV-faint evolved galaxies, 
which arc underrepresented in high redshift samples se- 
lected in optical surveys and are too faint for optical spec- 



troscopy. Second, we used the rest-frame optical spec- 
trum to study stellar populations and dust properties of 
galaxies, by modeling the spectrum and using spectral di- 
agnostics. Spectral modeling is very effective in narrow- 
ing the range of possible solutions for age, star formation 
timescale and dust content. Including NIR spectra in 
photometric modeling leads to more tightly constrained 
properties of stellar populations in high redshift galax- 
ies than modeling solely the optical-to-NIR photometry. 
And third, the break allows us to constrain the contri- 
bution from a possible AGN. In summary, studying the 
rest-frame optical shape appears to be crucial for under- 
standing the red z > 2 galaxy population, and provides 
information which is difficult - and in some cases impos- 
sible - to obtain by other means. 

The spectral modeling presented here for the 
three galaxies confirms the large variation among 
the SEP properties a nd S F Hs of DRG s (e.g. 
Forster Schreiber et all 120041: lLabbe et all 12003 
Webb et al.l 120061) . as our objects range from a dusty 
starburst with a small break to an apparently evolved 
galaxy with a strong 4000 A break and no detected Ha 
emission. For the galaxy without emission lines, we 
derived the redshift from the spectral continuum shape. 
The stellar masses of the two most evolved galaxies (age 
~ 1.3-1.4 Gyr) are about 4 x lO n M , and thus they 
are among the most massive galaxies yet identified at 
these redshifts. Comparison to low redshift galaxies 
suggest that they probably evolve into galaxies with 
properties comparable to those of the most massive 
and oldest galaxies in the low-redshift Universe. The 
younger starburst galaxy (0.3 Gyr) has a stellar mass 
of 1 x 10 11 Mq. The strong breaks in the two most 
evolved galaxies imply that stellar light is the dominant 
contributor to the optical emission. For the younger 
galaxy a significant contribution from an AGN cannot 
be ruled out. 

The next step to increase our understanding of high- 
rcdshift galaxy populations is to extend our NIR spec- 
troscopic sample. The present sample does not allow 
any statistical analysis, and is mainly an illustration of 
the capabilities NIR spectroscopy has to offer. A larger 
sample will tell us more about the ratio of evolved ver- 
sus star-forming z > 2 galaxies, and the contribution of 
AGNs. 



This research would have been impossible without the 
generosity and flexibility of the staff of the Gemini Obser- 
vatory. We are grateful to the referee for very useful com- 
ments, that greatly improved the paper. We thank Dawn 
Erb for providing the photometry and model parame- 
ters of the BX/MD sample. This research was supported 
by grants from the Netherlands Foundation for Research 
(NWO), and the Leids Kerkhoven-Bosscha Fonds. Sup- 
port from National Science Foundation grant NSF CA- 
REER AST-0449678 is gratefully acknowledged. DM is 
supported by NASA LTSA NNG04GE12G. EG is sup- 
ported by the National Science Foundation under Grant 
No. AST-0201667, an NSF Astronomy and Astrophysics 
Postdoctoral Fellowship (AAPF). ST acknowledges the 
support of the Danish Natural Research Council. The 
Keck Observatory was made possible by the generous 
financial support of the W.M. Keck Foundation. The 



12 



Kriek et al. 



authors wish to recognize and acknowledge the very sig- 
nificant cultural role and reverence that the summit of 
Mauna Kea has always had within the indigenous Hawai- 



ian community. We are most fortunate to have the op- 
portunity to conduct observations from this mountain. 



REFERENCES 



Adelberger, K.L., Steidel, C.C., Pettini, M., Shapley, A.E., Reddy, 

N.A., & Erb, D.K. 2005, ApJ, 620, L75 
Balogh, M.L., Morris, S.L., Yee, H.K.C., Carlberg, R.G., & 

Ellingson, E. 1999, ApJ, 527, 54 
Beers, T.C., Flynn, K., & Gebhardt, K. 1990, A J, 100, 32 
Bruzual A.G., 1983, ApJ, 273, 105 
Bruzual, G. & Chariot, S. 2003, MNRAS, 344, 1000 
Calzetti, D., Armus, L., Bohlin, R.C., Kinney, A.L., Koornheef, J., 

& Storchi-Bergmann, T. 2000, ApJ, 533, 682 
Daddi, E. et al. 2004b, ApJ, 600, L127 
Erb, D.K., et al. 2003, ApJ, 591, 101 

Forster Schreiber, N.M. et al. 2004, ApJ, 616, 40 

Forster Schreiber, N.M. et al. 2006, AJ, in press ( astro-ph/0510186 I 
Francis, P.J., Hewett, P.C., Foltz, C.B., Chaffee, F.H., Weymann, 

R.J., & Morris, S.L. 1991, ApJ, 373, 465. 
Franx, M. et al. 2003, ApJ, 587, L79 
Gawiser, E. et al. 2005, ApJS, 162, 1 
Giavalisco, M. et al. 2004 ApJ, 600, L93 
Hamilton, D. 1985, ApJ, 297, 371 
Kauffmann, G. et al. 2003a, MNRAS, 341, 33 
Kauffmann, G. et al. 2003b, MNRAS, 341, 54 
Kennicutt, R.C. 1998, ARA&A, 36, 189 
Kroupa, P. 2001, MNRAS, 322, 231 
Labbe, I. et al. 2005, ApJ, 624, L81 

Madau, P., Ferguson, H.C., Dickinson, M.E., Giavalisco, M., 

Steidel, C.C., & Fruchter, A. 1996, MNRAS, 283, 1388 
McLean, I. S., et al. 1998, Proc. SPIE, 3354, 566 
Papovich, C, Dickinson, M., & Ferguson, H.C. 2001, ApJ, 559, 620 
Papovich, C, et al. 2006, ApJ, in press ( astro-ph/0511289 I 
Pettini, M., Shapley, A.E., Steidel, C.C, Uuby, J.-U., Dickinson, 
M., Moorwood, A.F.M., Adelberger, K.L., & Giavalisco, M. 2001, 
ApJ, 554, 981 

Reddy, N.A., Erb, D.K., Steidel, C.C, Shapley, A.E., Adelberger, 
K.L., & Pettini, M. 2005, ApJ, 633, 748 



Rubin, K.H.R., van Dokkum, P.G., Coppi, P., Johnson, O., Forster 
Schreiber, N.M., Franx, M., & van der Werf, P. 2004, ApJ, 613, 
L5 

Rudnick, G., et al. 2001, AJ, 122, 2205 
Rudnick, G., et al. 2003, ApJ, 599, 847 
Salpeter, E.E. 1955, ApJ, 121, 161 
Sawicki, M., & Yee, H.K.C 1998, AJ, 115, 1329. 
Shapley, A.E, et al. 2005, ApJ, 626, 698 

Shapley, A.E., Erb, D.K., Pettini, M., Steidel, C.C, & Adelberger, 

K.L. 2004 ApJ, 612, 108 
Shapley, A.E., Steidel, C.C, Adelberger, K.L., Dickinson, M., 

Giavalisco, M., & Pettini, M. 2001, ApJ, 562, 95 
Simpson, C, Dunlop, J.S., Eales, S.A., Ivison, R.J., Scott, S.E., 

Lilly, S.J., & Webb, T.M.A. 2004, MNRAS, 353, 179 
Steidel, C.C, Giavalisco, M., Pettini, M., Dickinson, M., & 

Adelberger, K.L. 1996a, ApJ, 462, L17 
Steidel, C.C, Giavalisco, M., Dickinson, M., & Adelberger, K.L. 

1996b, AJ, 112, 352 
Steidel, C.C, Adelberger, K.L., Shapley, A.E., Pettini, M., 

Dickinson, M., & Giavalisco, M. 2003, ApJ, 592, 728 
Swinbank, A. M., Smail, I., Chapman, S., Blain, A.W., Ivison, R.J., 

&; Keel, W.C 2004, ApJ, 617, 64 
Tremonti, C.A. et al. 2004, ApJ, 613, 898 
Vanzella, E. et al. 2005, A&A, 434, 53 
van Dokkum, P.G. 2001, PASP, 113, 1420 
van Dokkum, P.G. et al. 2003, ApJ, 587, L83 
van Dokkum, P.G. et al. 2004, ApJ, 611, 703 

van Dokkum, P.G., Kriek, M., Rodgers, B., Franx, M., & Puxlcy, 

P. 2005, ApJ, 622, L13 
van Dokkum, P.G., et al. 2006, ApJ, 638, L59 
Webb, T.M.A., et al. 2006, ApJ, 636, L17 



